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We present a simple algorithm to implement the generalized derivative expansion introduced 
previously by L-H. Chan, and apply it to the calculation of the one-loop mass correction to the 
classical soliton mass in the 1+1 dimensional Jacobi model. We then show how this derivative 
expansion approach implies that the total (bosonic plus fermionic) mass correction in an N=l 
supersymmetric soliton model is determined solely by the asymptotic values (and derivatives) of 
the fermionic background potential. For a static soliton the total mass correction is — m/(2-7r), in 
agreement with recent analyses using phase-shift methods. 
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The calculation of quantum corrections to classical soliton masses is a key ingredient in the semiclassical approach 
to quantum field theory |l],|| . This question has been re-addressed recently [||4| for 1 + 1 dimensions using topological 
boundary conditions and phase-shift methods. In this Letter I compute the quantum mass correction for solitons 
in the recently introduced Jacobi model ||, using a simple algorithm based on the derivative expansion and I 

apply this derivative expansion approach to the total mass correction in an N — 1 SUSY model in 1 + 1 dimensions. 

Consider the following Lagrangian for a real scalar field <j> in 1+1 dimensions: 



V{4>) 



(i) 



We shall consider the examples: 



V{<t>) = { 



8A 



(VX ( 

V m 



l 



cos 



; <j) model 

; SineGordon model 
; Jacobi model 



(2) 



The </> 4 and Sine-Gordon models are text-book cases 0] , while the Jacobi potential has recently been studied in the 
context of an exactly solvable model of quantum mechanical instantons Since quantum mechanical instantons 
have the same functional form as solitons in 1+1 dimensions, some of the results of j5| carry over here. The function 
sn(z|^) is one of the Jacobi elliptic functions ||, and < v < 1 is the (real) elliptic parameter. This parameter v 
controls the shape and period of the Jacobi potential, as illustrated in Figure 1. The Jacobi model is a deformation 
of the Sine-Gordon model, reducing smoothly to Sine-Gordon as v — ► [note that sn(^|y = 0) = sin(§)]. 

As is well known jj]], the mass m and coupling A may be scaled out, so that the semiclassical loop expansion is an 
expansion in powers of -^4. With this understood, we can set m = A = 1 in our one-loop calculations. 

These models have classical static soliton/antisoliton solutions 4> c (x) satisfying 



<f/ e (x) = ±^/2V{4>c{x)) 
For the potentials in ^ the classical solitons are 

)tanh(|) ;</> 4 model 
4arctan(e a: ) ; SineGordon model 
2K(v) + 2 sn _1 (tanh(a;)) ; Jacobi model 



(3) 



(4) 



Here K(v) = J*' dO/y/l - t/sin 2 6», is the elliptic quarter period f§. Note that K(0) = tt/2, and K(u) - ^log(^) 
as v — ► 1. 
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FIG. 1. The Jacobi potential V{4>) = 2sn 2 (||^). The solid curve is for v = 0.99, while the dashed curve is for u = 0, which 
is just the Sine-Gordon case. 

The classical mass of the soliton is given, in units of ^r-, by 

f f ; 4 model 

Af classical = / dty/Wtf) =\ 8 : SineGordon model (5) 
J ™ ia l^ulog(^§J ; Jacobi model 

where the integration is between the two neighboring minima of V(<f>) between which the soliton interpolates. 
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FIG. 2. The bosonic fluctuation potential V" (<j) c (x)) for the Jacobi model. The solid curve is for v — 0, which is the same 
as the Sine-Gordon case. The dashed curve is for v = 0.6, and the dotted curve is for v — 0.99. When v 7^ 0, this fluctuation 
potential is not of the solvable Poschl- Teller form. 



In units of m, the one-loop quantum correction to this classical soliton mass is 

(6) 



M = - / — tr log 



2tt 



This expression must be suitably renormalized, but once this is done M is a pure number. Here V" (<p c (x)) means 
-£jpV (4> c (x)) . For the cases listed above, the bosonic fluctuation potential V"((p c (x)) is 

fl-§sech 2 (§) ;0 4 model 
V"(<f> (x)) = < l-2sech 2 (x) ; SineGordon model (7) 

Id--) ( - 1+ V^ + r <w ) ; Jacobi model 

For the </> 4 and Sine-Gordon cases, the fluctuation potentials in ([?]) are of the exactly solvable Poschl- Teller form, 
which means that the quantum mass correction in (pi) can be computed exactly. The fluctuation potential for the 
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Jacobi model is shown in Figure 2; its spectral properties are not known exactly, and so we need an approximate 
method to compute the mass correction (ph. Here we use the derivative expansion ^,0]. 

To implement the derivative expansion we first note that V" (<fi c (x)) — > 1 as x — > ±oo, for the fluctuation potentials 
appearing in (^). Moreover, this background value of 1 (equal to m 2 with mass scales reinstated) is approached 
exponentially fast. So it is more natural to expand in terms of the difference 

W{x) = V"{<p c {x)) - 1 (8) 

and its x derivatives. This leads to a generalized derivative expansion Q in which (^) is rewritten and expanded as 

M [W] = - f ^c 2 tr ( 1 ) 

= [°° ¥-u*ziA(z,W\ (9) 

where z = 1/(uj 2 + 1). In the derivative expansion approach ^M, the expression @ is renormalized by dropping a 
divergent term that is independent of W (this cancels in taking the difference with the W = case) , and a divergent 
term linear in W (this corresponds to making the bosonic tadpole graph vanish, as in Then A(z 7 W] is defined 

by the expansion 

-. oo 

A(z,W}=~-zY / z n anlW] (10) 

71=0 

where the expansion coefficients a„[VT] are functionals of W{x), the first few of which are 0: 
a [W] = - / W 2 dx 

° J-oo 
c poo 

ai[W} = -^- {2W 3 + (W / ) 2 } dx 

^ Z J-oo 

7 f°° 

a 2 [W] = — / {5W 4 + 1QW(W') 2 + (W") 2 } dx 
128 J_ oc 

q r°° 

a 3 [W} = -— {UW 5 + 70W 2 {W') 2 + UW(W") 2 + (W'") 2 } dx (11) 
512 J-oo 

Given the expansion ( |l0| ) for A(z, W], the lo integrals in (^|) can be performed to yield a simple formula for the soliton 
mass correction: 

M [W] = - J= Y X {n + a n [W] (12) 

Thus, the derivative expansion computation of the soliton mass correction reduces to the calculation of the a„[W] in 
( |TT| ) for the given W(x) in (||). We stress again that each of these a„[V7] is a pure number |p^| . 

The functionals a^fW] can be obtained by a direct expansion of the resolvent in (||), as in M. However, it is 
easier computationally to note that these functionals have simple recurrence relations (these functionals also arise in 
asymptotic and WKB expansions of discriminants [p|JIo|, and in zeta function analysis (ll)). In fact, 



a n [W] = (-l) n+1 (2n + 3) / r 2n+4 (x)dx (13) 

J — oo 

where the functions r n (x) are defined by the following recurrence relation: 
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r Q (x) = l , r x (a;) = , r 2 (x) = ~W(x) 



i 1 

J=l 

This algorithm for generating the r2 n +4:(x) functions is simple to implement using, for example, Mathematica or 
Maple. Furthermore, all the integrals in ( |13| ) can be performed analytically for each of the three soliton models listed 
in®. 

So far, this algorithm could be applied to any suitable potential in (^), with W(x) — > sufficiently quickly 

as x — > ±oo. If the background is that of a soliton <j> c {x) satisfying (JsJ) , so that W(x) is given by (||), it is more 
convenient to change variables in the integration in ( |l3| ) from x to <j) c , using (Jfy. For example: 

a [W] = *[ -^W* 
8 J min 

a2 [ W ] = JL f J£ \5W 4 + 20WV(V^) 2 + (V'V^ + 2VV&) 2 } 

o 3 [W] = _A / W \lAW 5 + UOW 2 V(V^) 2 + UW(V'V^ + 2VV^f 
512 J min y/2V((j)) I- 

+2t/(t/"y (3) + 3VV (4) + 2VV^ (5) ) 2 } (15) 

These modified expressions for a n [W] are straightforward to generate, either by converting the expressions in (^), or 
by a new recurrence relation based on (fl4"|). In (15), both V and W = V" — 1 are treated as functions of <j>, and the 
integration is between the neighboring minima of V(cj>) between which the soliton interpolates, just as in the standard 
calculation (|J) of the classical mass. The advan tag es of this form (JlJ) are : (i) the only input is the original potential 
V{4>)\ (ii) the integrals are easier than those in (|ll|); (iii) it organizes the expansion ( |l2| ) as a derivative expansion, as 
each a n in ( |l5| ) involves only terms with (2n + 4) derivatives of V with respect to 4>. 

The derivative expansion calculation of the soliton mass correction ([l2]) can now be programmed as a simple 
algorithm: generate the expressions in ( |ll"| ) or (|l5|) using recursion relations, do the integrals, and then insert the a n 
into (Q). We first test this algorithm on the well-known solvable cases. For the Sine-Gordon model one finds a^p = 2 
for all n, confirming a result in 0. Then, using our formula (|l2|), we immediately obtain the standard result [l]j|] 

,.sr 1 r(n+l) 1 

M = > — — = — fl6) 

For the </> 4 model, one finds af^ = 2 + (l/4)™ +1 for all n, again confirming a result in 0. Then, from ( [12] ) we find the 
standard result |l|,|[ 



1 _ _±_ T(n+l) = J_ _3_ , m 

7T 320F ^ 4™r(n + 5/2) 4\/3 2vr 1 ' 

One can speed up the convergence of the derivative expansion by separating out the contribution from the zero mode 
of the fluctuation potential [j7|. This is achieved by writing A(z, W] = (A(z,W] + ^) — in (Q) and (|lO|), which 
separates out a — — term, so that the mass correction (|l2| ) becomes 

M[W] = -i - -1= f] ^"V/l («nW - 2 ) ( 18 ) 
7T s-v/tt ' 1 (n + 5/2 

v n— 

In this form (Il8|), already the leading (n = 0) derivative expansion term gives the exact answer for the Sine-Gordon 
case. For the cjr model, the leading term gives 10% accuracy, the first correction gives better than 2% accuracy, two 
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corrections give 0.3% accuracy, and three corrections give 0.05% accuracy. This is impressive accuracy, especially 
considering ao, ■ . ■ , as can easily be calculated analytically. 

For the Jacobi model there does not appear to be a simple closed formula for the a n . This is presumably a reflection 
of the fact that in this case, unlike in the c/> 4 and Sine-Gordon cases, V"(<j> c (x)) is not an exactly solvable potential. 
Nevertheless, we can compute various orders a n of the derivative expansion expression with ease: 



^ Lame 



Lame 



128i^ 3 / 2 

1 



2V^(5 - 38v + hv 2 ) - (5 + Zv + 3v 2 + 5v 3 ) log 

. .„ [-2VR1785 + 520^ - 11746^ 2 + 520j/ 3 + 1785i/) 
12288^ 5 / 2 ' 



( 1+ 


^1 







-15(119 - 5i> + Uv 2 + Uis 3 - 5v 4 + ll9v 5 ) log 



1 - V" 



1 



[-2^(1323945 - 683970^ + 146839^ 2 - 3800572^ 3 + 146839i/ 4 - 683970^ 5 + 1323945z/ 6 ) 



3932160^ 7 /2 

+ 105(12609 - 10717^ + 2449^ 2 + 779v z + 779v 4 + 2449t/ 5 - 10717^ 6 + 12609;/) log 

- rg70no L n 9/2 [2^(258321525- 355550860^+ 143327940J/ 2 - 681332zy 3 

-194557554t/ 4 - 681332i/ 5 + 143327940;/ - 355550860z/ + 258321525z/) 
-105(2460205 - 4206267^ + 2275076i/ 2 - 275052i/ 3 - 24586i/ 4 - 24586i/ 5 

-275052i/ 6 + 2275076i/ 7 - 4206267:/ + 2460205i/ 9 ) log 



(19) 



Notice the appearance of the classical soliton mass log ^ jrr^= J in these expressions for the one- loop quantum 
correction. 

Figure 3 shows the successive derivative expansion approximations to the quantum mass correction, as a function 
of the elliptic parameter v. The expansion appears to be convergent for each v. Furthermore, a nontrivial test of the 
results in ([H]) is that they must reduce to the Sine-Gordon results as v — > 0. Indeed, from ( |l9|) we verify that for 
each n 



as v 







(20) 



On the other hand, as v — > 1 the quantum mass correction diverges, just as does the classical Jacobi soliton mass in 
(§). This singular behavior is due to the fact that as v — > 1 the period of the Jacobi potential diverges logarithmically, 
so neighboring vacua become infinitely separated, and the asymptotic limits of the classical soliton (f> c (x) in (Q) diverge. 
We now consider the N — 1 supersymmetric extension of the bosonic soliton model ([IJ) . Consider 



£susy = \{d^){d»<t>) ~ \u 2 {^) + - \u\4>)^ 



(21) 



where tp is a Major ana fermion and 

u{<t>) = 

The bosonic potential is V(4>) = 2 



-2 Bin(|; 



1) 



-2sn 



(fl") 



; </> 4 model 
; SineGordon model 
; Jacobi model 



(22) 



-U 2 (<p), solitons are given by — 



dx 



v F {4>) = u'{4>) 



-U((j> c (x)), and the boson-fermion coupling is 

(23) 
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FIG. 3. The successive derivative expansion approximations for the quantum mass correction in (jH^) for the Jacobi model, 
using the coefficients a n in ([H)]), plotted as a function of the elliptic modular parameter v. The dashed curve is the leading 
(n = 0) contribution, the dotted curve includes the first two terms, while the solid curve includes the first three terms, and is 
indistinguishable on this scale from that with the first four terms. The solid horizontal line is the Sine-Gordon mass correction 
M SG = — — , to which the Jacobi model reduces when u = 0. 

The fermionic contribution to the one-loop soliton mass correction is 

Mp = i tr log[^- V F (c/> c {x))] 

= -- {M[W+\ + M[W-]} (24) 

where M[VF] is the bosonic mass correction in (||) and 

W±(x) = [Vp^ix))] 2 ± ^V F (4> c (x)) - 1 (25) 

ax 

But with the SUSY coupling (p3|), W_(x) is equal to the bosonic potential W(x): 
Thus the total one-loop mass correction, including both bosonic and fermionic contributions, is 



W^(x) = W{x) = —V(cp c (x))-l (26) 



M tota i = i {M[WJ\ - M[W+]} (27) 

That is, the net mass correction is (half) the difference of two bosonic mass corrections, for W^{x) respectively. 

Now notice the remarkable fact that if we compute this difference using the derivative expansion (|l2|), then term 
by term each of the functionals a„[W_] — a^W+J obtained from ( [TT| ) is an integral of a total derivative. For example, 



a [W-] 


-a [W+ 


ai [wq 


-ai\W+ 


a 2 [W-} 


-a 2 [W+ 




-a 3 [W+ 



m /I f Tx { 7Vf ~ 7Vf + T Vf Vf + 7{Vf Vf){Vf)2 ~ s W*> a " Yo Vf ^ )2 } dx 

-ikFJ^i {~ Wf + 12Vf " t v * + y^ 7 _ Vf 18Vf{1 ^ m)2 _ WFiVF)i 



~(i - v 2 )(v F fv F < + ™v F (i - v'w^f + g(u;') 3 - w - ^v F (v F y } <i>- 



(28) 



G 



In (f28|), Vp means -^Vf(4> c (x)). It is straightforward to generate these expressions to a desired order using Mathe- 
matica or Maple. 

Due to this dramatic simplification, the differences a n [W_] — a n [W_|_] are determined solely by the asymptotic values 
(at x = ±00) of Vf{x) and its x derivatives. But for a solitonic background, Vf(4>c(x)) approaches ±1 exponentially 
fast, so 



Vf(4> c (x)) = ±1 at x = ±00 

A'- 

■ F 

For example, these conditions are clearly satisfied for the (f> 4 , Sine-Gordon and Jacobi models, 



V F n) {<Pc{x)) = at x = ±oo , n>l (29) 



!tanh(|) ;0 4 model 
tanh(a;) ; SineGordon model (30) 
(i-i/)tanh(a:) Jacobi model 
1— vtanh (x) ' 



Inserting the asymptotic conditions (|29j) into (|28|), the derivatives of Vp do not contribute, and we discover that the 
differences are independent of n: 

On[W-] - a n [W+] = 2 , for all n (31) 

Thus, combining the derivative expansion formula ( |l2"| ) with the expression ( p7[ ) for the total mass correction, we find 
(in units of m) 

_ 1 ^ r(n + l) _ 1 
M total - ^ r(n + 5/2) - (32) 

v n— 

This is in agreement with recent analyses ||,|J using topological boundary conditions and phase shift methods. We 
note that in the derivative expansion approach, this result is quite general, relying only on the asymptotic values of 
Vp(4> c (x)) and its x derivatives. 
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